Numerical study of Fermi-Pasta-Ulam recurrence for water waves over finite depth 
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Highly accurate direct numerical simulations have been performed for two-dimensional free-surface 
potential flows of an ideal incompressible fluid over a constant depth h, in the gravity field g. In 
each numerical experiment, at t = the free surface profile was in the form y = Ao cos(27rx/L), 
and the velocity field v = 0. The computations demonstrate the phenomenon of Fermi-Pasta-Ulam 
(FPU) recurrence takes place in such systems for moderate initial wave amplitudes Aq < Q.12h 
and spatial periods at least L < 120h. The time of recurrence Tppu is well fitted by the formula 
TFPv{g/hy/^ « 0.16{L/hf{h/Aoy/^ 

PACS numbers: 47.35.-1, 47.15.K-, 47.10.-g 



Many nonlinear dispersive waves are known to exhibit 
the Fermi-Pasta-Ulam (FPU) recurrence, when a (finite- 
size) system nearly repeats its initial state after some pe- 
riod of evolution. For the first time, this phenomenon 
was observed in the famous numerical experiment [l| 
with one-dimensional (ID) lattices of nonlinear oscilla- 
tors. The impact of that observation on the subsequent 
development of nonlinear science was very deep. In par- 
ticular, it resulted in the discovery of solitons by Zabusky 
and Kruskal [2]]. Since then, the FPU recurrence and re- 
lated phenomena were studied in many physical contexts 
(see, e.g., [3l-[T3|. and references therein). 

In general, the FPU recurrence takes place in a system 
if its dynamics is nearly integrable. Two different reasons 
can result in such approximate integrability. In the first 
case, just a few (2-3) collective degrees of freedom are 
effectively involved into the evolution due to a small sys- 
tem size and low level of nonlinearity, and that restricted 
dynamics is integrable (for example, such variant of FPU 
recurrence takes place for 3D deep-water waves [tI. ITi|). 
In the second possible case, the evolution is governed by 
equations of motion which are close to some completely 
integrable system The first case does not require 
any additional symmetry and therefore is more common, 
while the second case is more interesting, since it allows 
for a large number of degrees of freedom to be involved 
into the process. 

As the theory of water waves is concerned, there are 
three the most popular integrable models for three es- 
sentially different dynamic regimes: (i) the nonlinear 
Schroedinger equation (NLSE) approximates an envelope 
of deep-water waves, (ii) the Korteweg-de-Vries (KdV) 
equation describes weakly nonlinear dispersive unidirec- 
tional shallow- water waves, and (iii) the Boussinesq equa- 
tions approximately describe bidirectional shallow-water 
waves (concerning integrability of the Boussinesq equa- 
tions, see While the FPU recurrence was in- 
vestigated in detail for the regimes (i) IHI and (ii) 
, it was not studied for the regime (iii) , which case in- 
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FIG. 1: Three examples of FPU recurrence in water wave 
systems: the free surface profiles at f = and ait — Tfpu(/i = 
1 m; L; Aq = 0.12/i) are very close for L/h — 40 and for L/h = 
60; for L/h = 80 the recurrence is not so perfect but still 
apparent. The curves with the highest elevation correspond 
to the moments when a "soliton" collides with the wall. 



eludes water waves propagating in a closed flume. The 
present work is intended to fill this gap by means of direct 
numerical simulations. 

It is important that nonlinear dynamics of free water 
surface is not exactly integrable (see, e.g., and ref- 

erences therein) , and therefore it was not obvious a priori 
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FIG. 2: The maximum and minimum elevations of the free 
boundary versus time for L/h = 40 and different initial am- 
plitudes. The peaks of the dependences at approximately reg- 
ular time intervals correspond to collisions of "solitons" with 
the "walls". 



FIG. 3: FPU recurrence for Ao/h = 0.12 and different L. 
With increasing L, the dynamics between the recurrence cases 
becomes more and more complicated: additional minima in 
the envelope of Ymax appear, which correspond to significant 
excitation of the second, or the third, or a higher Fourier mode 
in nearly standing-wave regime. 



in what parameter region the FPU recurrence can be ob- 
served. To clarify this question, here the highly accurate 
numerical method based on exact equations of motion 
for two-dimensional (2D) free-surface potential flows of 
an ideal incompressible fluid was employed (20I . [2l| . in 
the most simple variant when the bottom is horizontal 
at a constant depth h. In each numerical experiment, 
at i = the free surface profile was taken in the form 
y ^ Aq cos[2'kx/ L) (the gravity acceleration g is directed 
against y axis), and the velocity field was zero every- 
where. Thus, this configuration corresponds to waves in 
a basin with vertical walls at x = and x — L/2. 

Since the exact equations of motion are written in 
terms of so called conformal variables, and the surface 
shape is given in a parametric form X + lY = —ih + 
(1 + iR)p{d,t), with R being a linear integral operator 
[13, [in, it was a nontrivial task how to determine the real 
function p('i?, 0) corresponding to a given initial surface 
profile y = rio{x). This technical problem was solved by 
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FIG. 4: The time of recurrence, excluding tire cases 
when {Ao/hy/'^{L/h) < 1, is well fitted by the formula 
Tppvig/h)'^/'^ « C{L/hf{h/Aof^'^, with C = 0.148 + 
0.096(Ao//i). Roughly, C = 0.16 ± 0.01. 



an auxiliary numerical procedure, when in one of the two 
evolutionary equations (namely, in the Bernoulli equa- 
tion) the combination {ipt + gY} was temporarily re- 
placed with the combination {ipt + 9[Y — ^70(2;)] + T?/'}, 
where ^{'d,t) is the surface value of the velocity poten- 
tial, and r — const > is some artificial damping. With 
this modification, the free surface quickly evolved from 
y = to the required initial profile, a.nd after that we 
turned to the original equations [20, [21| . 

The parameters Aq and L were taken inside the re- 
gion 0.04 < A^/h < 0.24 and 20 < L/h < 120. The 
computations have demonstrated that the phenomenon 
of FPU recurrence indeed takes place for moderate ini- 
tial wave amplitudes Aq < 0A2h. After a few oscil- 
lations in a standing-wave regime, the system enters a 
regime with soliton-like coherent structures moving be- 
tween the "vertical walls" . The shapes of the "solitons" 
evolve with time, and after a period Tppu, which depends 
on h, L, and Aq, the system again enters the standing- 
wave regime and approximately repeats the initial state, 
as it is shown in Fig.l. With a realistic value h ^ 1 m, 
the period of recurrence is several minutes. 

In Fig. 2 and Fig. 3, the maximum and minimum eleva- 
tions of the free boundary versus time are shown for dif- 
ferent Aq and L. These dependences clearly indicate the 
presence of the FPU phenomenon. However, a "quality" 
of the recurrence is not uniform in the parametric region: 
it is better for smaller initial amplitudes and shorter spa- 
tial periods L. For large amplitudes Ao/h > 0.20, one 
can hardly recognize more than one recurrence case. In 
that strongly nonlinear regime, very high waves grow in 
the system, which can eventually produce sharp crests 
and break (not shown). 



The time of recurrence rapidly increases with the spa- 
tial period L, and behavior of the system between the 
recurrences becomes more complicated, since a larger 
number of solitons participates in the dynamics, Ng ~ 
{L/h){Af)/hY/'^ . In the course of evolution, there are 
some time intervals when a significant part of the energy 
is concentrated in the second, or in the third, or in a 
higher Fourier mode. 

The numerical results for Tppu for different L and Aq 
are summarized in Fig. 4, which shows that, excluding the 
cases when {L/h){Ao/hy^^ < 1 (under this condition, 
the length L becomes too small to contain a soliton with 
an amplitude about ^0)1 the time of recurrence is well 
fitted by a formula 



Tfpu (f ) 



1/2 



0.148 



0.0964^ 

h 



h 

A'o 



1/2 



(1) 

It should be noted here that in each case the quantity 
Tppu was determined rather roughly, just by looking at 
graphs like those shown in Fig. 2 and Fig. 3. Therefore 
the corresponding accuracy is rather low, despite the fact 
that the accuracy of the simulations is very high (so, the 
energy conservation was at least 7-8 decimal digits from 
the beginning of each computation to its end) . 

In order to observe the FPU recurrence in a real-world 
flume filled with water, the viscous friction near the bot- 
tom and near the side walls should be taken into account. 
Thus an additional condition arises, 7 • Tppu ^ 1, where 
the damping coefficient 7 can be estimated as 



(2) 



where ly is the kinematic viscosity. Since with fixed ratios 
Ao/h and L/hthe time of recurrence behaves as h^^^, the 
product 7 • Tppu behaves as h~^^'^, and therefore it can 
be indeed small with sufficiently large h. 

To conclude, in the present work the phenomenon of 
Fermi-Pasta-Ulam recurrence was observed in the numer- 
ical experiments modeling potential one-dimensional wa- 
ter waves in a closed basin. Based on the numerical re- 
sults, a simple fitting formula for the period of recurrence 
is suggested. 
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